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A model of the magnetic vortex in a thin disk of soft magnetic material is constructed. The model 
employs a piecewise merger of two of the most successful analytical models of the vortex state in 
a disk and is capable of describing the change in the characteristic warping of the magnetization 
distribution when the vortex core interacts with a magnetic pinning site. The constructed analytical 
model is compared to numerical simulations of ideal disks with and without pinning sites and 
is found to accurately predict the magnetization, vortex position, hysteretic transitions, and 2-D 
displacement of the vortex in the presence of pinning sites. 

PACS numbers: 75.30 Hx, 75.60.Ch, 75.75.Fk 



I. INTRODUCTION 

Interest in magnetic vortices^ in thin disks has grown 
dramatically over the past two decades, as these are fun- 
damental physical systems with direct applications to 
technology 3 . Topological structures such as vortices are 
stable, manipulable objects that show promise as logic el- 
ements or storage media in spintronics applications. The 
thin soft magnetic disk, the prototypical system contain- 
ing a vortex, has therefore become an extensively inves- 
tigated system. Properties studied include struc ture^, 
dynamical modea^^, annihilation™^ and creatiorPE^ISl. 
As each aspect of vortex physics is probed experimen- 
tally, and considered for technological applications, the- 
oretical understanding via simulation and modeling is 
also advanced. Modeling is particularly important in the 
case of the thin ferromagnetic disk as it presents a well- 
defined system amenable to description by an analytical 
approach. 

A recent focus of activity is the interaction of vortex 
cores or domain walls with film inhomogeneities. Phys- 
ically, geometric defects or magnetic impurities can in- 
crease or decrease the energetic cost for the topological 
structure^ creating preferential locations or altering the 
magnetization distribution. In the disk system, direct 
observations of vortex core pinning have been observed 
with Lorentz microscop}El while the effect on vortex gy- 
ration has been observed with time-resolved magneto- 
optical Kerr effect microscop}^^ and electronical!}!^. 
Incorporation of pinning potentials into existing analyt- 
ical models has permitted a qualitative description of 
the position of the vortex and its reduced displacement 
susceptibility 21 . This approach is insufficient for quan- 
titative applications. Recent work using nanomechan- 
ical torque magnetometry has provided direct observa- 
tion of the Barkhausen steps associated with jumps in 
core position 2 ^, necessitating the development of a model 
that permits a quantitative description of pinning ef- 
fects. Physically one expects two clear contributions to 
the magnetic susceptibility of a pinned vortex, one from 
the magnetization distribution and one from rigid dis- 
placement of the vortex. An analytical model designed to 



capture both contributions is presented here, and demon- 
strated to reproduce numerically simulated quasistatic 
vortex pinning behavior. 

The vortex state ansatz was first developed for the 
magnetic disk by Aharoni in lOOO 2 ^. Further work by 
Usov and Peschanny 4 determined an exchange optimized 
functional form of the core magnetization profile. Good 
agreement between this model and simulation 2 -^, as well 
as experimental observation 5 , was found. This work con- 
sidered the vortex ground state, at zero field. Subsequent 
workpHHl considered the displacement of the core with 
field and annihilation field using a model that rigidly 
translated the magnetization distribution developed by 
Usov and Peschany. This is known as the Rigid Vortex 
Model (RVM). Recent extensions of the RVM include 
higher order versions developed to describe the suscep- 
tibility of the displaced vortex analytically 10 . Concur- 
rently, more complex models approximating t he flex ing of 
the magnetization distribution were developecPzHH] This 
class of model is equivalent to the influence on the mag- 
netization distribution of a second vortex moving from 
infinity to the edge of the disk, and it is called the Two 
Vortex Model (TVM). Versions of this model have been 
applied to calculating the stability of the vortex state 
in a disk with moderate success^^S and, with greater 
success, to predict frequencies of dynamic modes^HI]] 

II. CONSTRUCTION OF THE MODELS 

Both the RVM and TVM link the computed model 
magnetization directly to the vortex core displacement 
in the disk using only a single parameter^. For each 
model the total energy is composed of the exchange, de- 
magnetization, and external field energies as summarized 
below. Magnetocrystalline anisotropy is neglected. 

A. The Rigid Vortex Model 

The RVM is derived by considering the vortex magneti- 
zation distribution ansatz^Uto be immutable and trans- 



2 




FIG. 1. (Color available in online version.) Schematics de- 
picting the evolution of the magnetization distribution as the 
vortex core is displaced. The arrows, as well as the color scale, 
indicate the in-plane direction of the magnetization. At top 
is the DVPM which has two regions: an outer shell described 
by the RVM and an inner shell described by the TVM. Both 
components make contributions to the magnetization and vor- 
tex displacement allowing the DVPM to capture more com- 
plex behavior. By contrast the RVM (lower left) employs a 
parameter that governs both vortex and magnetization dis- 
tribution displacement, while the TVM (lower right) has a 
single parameter that determines both magnetization distri- 
bution flexing and vortex displacement. 



to the second order RVIvP^with m (h) = b Q (h) = h//3. 



B. The Two Vortex Model 

The TVM is derived by setting a boundary condi- 
tion and computing the magnetization distribution that 
minimizes t he ex change energy for a given vortex core 
displacement! 27 1 35 1 (F igure ^p) . Here the no side charges 
version of the mode l 27 ! 29 ! is applied and the magnetiza- 
tion is held tangential to the disk edge. In the same form 
as used for the RVM, the total normalized energy may 
be written down for the TVM, 



tot 



a 



^M 2 V 2 a 



(2) 



where a/2 = Ar/R represents the normalized displace- 
ment, a = RFi(L, R)/L — R 2 Q /2R 2 incorporates the de- 
magnetization energy and exchange energy, and £ is a 
constant (~ 10/29). The function Fi(L,R) is an equiva- 
lent demagnetization factor describing the volume mag- 
netostatic charges resulting from flexing of the magneti- 
zation distribution and is approximated as k(L/R) 2 with 
k = 0.08827. As before, minimization with respect to a 
allows computation of a Q (h) = ^h/a and magnetization 

™>o(h) = &o(h). 



lating that distribution relative to the physical boundary 
of the disk (Figure [l]). The normalized total energy for 
the 3rd order RVM 1Tr of a disk with a radius R and thick- 
ness L as a function of the reduced field h = H/Ms is 
given as 

where b = Ar/R is the normalized core displace- 
ment, f3 = F C (L,R) — R 2 /R 2 is a constant describing 
the demagnetization energy and exchange energy with 
F C (L, R) representing the susceptibility-corrected demag- 
netization factor com puted for t he uniformly magnetized 
diskP^, and R Q = ^2A/ \i Q M^ is the exchange length. 
The energy is normalized by /i G M|V where V is the 
disk volume. It should be noted that the incorpora- 
tion of the susceptibility correction to the demagneti- 
zation factor, the so called /i* correction^, is critical to 
the success of the model in application to a disk com- 
posed of a soft magnetic material (Figure [2] a, inset), 
but is often neglected. Solving the model permits calcu- 
lation of the vortex displacement as a function of field, 
b (h) = (-4/3 + 2^4/32 + 6ft 2 )/3ft, and consequently the 
magnetization m (h) = M (h)/Ms = b Q — b^/8. Re- 
moving the third order term in the magnetization, and 
consequently the external field energy, reduces the model 



C. The Deformable Vortex Pinning Model 

In the RVM and TVM models, the direct dependence 
of the magnetization on the vortex displacement param- 
eter means that, even in the presence of a local imperfec- 
tion in an otherwise perfect disk , the core position and 
overall disk magnetization can not evolve independently. 
This implies, and numerical simulation confirms, a clear 
qualitative shortcoming of these models in the form of the 
flexing of the magnetization distribution evolution dur- 
ing pinning (See Supplementary Movies Ml and M2). A 
composite model, the Deformable Vortex Pinning Model 
(DVPM), can be constructed to correct this problem, in- 
corporating both the RVM and TVM with two parame- 
ters. 

The DVPM can be derived by dividing the disk into 
two regions, an outer annulus described by the RVM, 
and an inner disk of radius R\ described by the TVM 
with no side charges. The tangential boundary con- 
dition of the TVM maintains a piecewise continuity of 
the magnetization distribution between the two models. 
As field is applied, the center of the inner disk shifts 
along with the core position computed from the outer 
annulus. Concurrently, the computed core position in 
the inner region also deflects, deforming the central re- 
gion magnetization. In this way, there are two contri- 
butions to vortex position as well as magnetization dis- 
tribution. To self-consistently solve the two sections of 
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the model, it is necessary to consider the energetic com- 
ponents of each (Figure [2| . Neglecting the susceptibil- 
ity correction, derivation of the inner TVM energetic 
contributions would be independent of the outer annu- 
lus, and since the RVM contributions change only from 
the shift of the magnetization distribution relative to the 
boundary, the two models would be independent up to 
the point where the central region edge reaches the disk 
boundary. However, the susceptibility correction implic- 
itly couples the two sections through the demagnetiza- 
tion energy. In the RVM, the susceptibility correction 
uses material susceptibility to estimate the net decrease 
of the demagnetization energy due to the reduction of 
side charges, and the introduction of volume charges pa- 
rameterizing the real world magnetization flexing. But 
the TVM is already flexible, and therefore including the 
susceptibility correction for the shell, in addition to a 
central TVM region reduces the demagnetization cost 
of core deflection too much. This can be compensated 
for by increasing the energetic cost of deflection for the 
rigid shell a corresponding amount. This is equivalent to 
reducing the magnitude of the susceptibility correction. 
A linear interpolation between the corrected (F c ) and 
uncorrected (F nc ) demagnetization factors of the form 
F(L, R, R t ) = (1 - R 1 /R)F C (L, R) + (R 1 /R)F nc (L, R) is 
introduced. The equivalent change of the RVM shell de- 
magnetization factor necessary to rebalance the energetic 
cost of vortex displacement is (Ri/R) 2 (2Ri/L)Fi(Ri, L). 
Comparing each, the deviation does not exceed 10% over 
the range R\/R =0 to 1 for a linear interpolation, cor- 
roborating this interpretation^ . 

Using the interpolated demagnetization factor, the en- 
ergy of the combined piecewise model may be written, 



^M 2 V ~ 2 



F«2 



h(b f 



12 



iho!\ (3) 



where = F{L,R,R{) - R 2 Q /R 2 and a' = 

RiF x {L, R^/L - R 2 /2R 2 . Here b' is the normalized dis- 
placement of the outer RVM shell, and a' /2 = Ari/Ri 
is the central TVM core displacement normalized to R\. 
The factor 7 = R 2 /R 2 scales the energy contributions ac- 
cordingly. All other symbols remain as before. The total 
core displacement is b' + R\o! j2R with the same expres- 
sions for b' (h) and a' (h) as before, but with j3' and a' 
replacing j3 and a respectively. The corresponding mag- 
netization is m (h) = b' {h) - 6 / (/i) 3 /8 + 7^ (ft). For 
small displacements, the third order terms for the RVM 
shell may be dropped to make a simplified version of the 
model applicable for small displacements. For an ideal 
disk, the choice of R\ is not critical, the susceptibility 
correction interpolation mitigates variation with R\ for 
a wide range of disk aspect ratios. A value R\ = R/2 is 
used below. 
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FIG. 2. (Color available in online version.) (a) The com- 
puted m-h curves from the four models discussed in this 
paper are compared against a simulation of a 1 /im diam- 
eter, 30 nm thick disk with M s =800kA/m. (b) The com- 
puted vortex displacement as a function of field for each 
model is compared against the simulation. The legend in- 
set in panel (b) applies to both (a) and (b). Only the RVM 
with a susceptibility-corrected demagnetization factor and the 
DVPM describe both position and magnetization accurately 
for displacement <R/2. Inset in (a) is a comparison of the 
interpolation (solid red) used for the susceptibility correction 
applied in the DVPM against the effective demagnetization 
factor required to balance the displacement contributions of 
the RVM and TVM regions in the DVPM. This correction 
makes the computation insensitive with respect to variation 
of the size of the TVM region (radius R x ). All DVPM calcu- 
lations in panels a-d use R\ =R/2. (c) The computed initial 
susceptibility for all models is compared against simulation 
for disks of varying size. Simulations denoted with triangu- 
lar point use a disk diameter of 500 nm, circular points use 
1 /xm, and triangles use 3.6 /im. All simulations use Ms — 
800kA/m except the 3.6 /im which use M s =715kA/m. (d) 
Using the same simulation parameters the mode frequency of 
the gyrotrpic mode was computed. The legend in (c) applies 
to panel (d) as well. Only the DVPM agrees well with both 
the initial susceptibility and the gyrotropic frequency. The 
corrected RVM correctly computes the susceptibility while 
the TVM closely estimates the gyrotropic frequency. For (c) 
and (d) calculations were performed holding R =500 nm with 
variable thickness. 



III. APPLICATION TO THE IDEAL DISK 

From Equations [I] [2j and [3] the ideal disk behaviour of 
each model may be computed and compared to Landau- 
Lifshitz-Gilbert micromagnetic simulation^ . To mimic 
quasistatic behaviour, time integration with a damping 
factor of 1.0 was used. All simulations were performed 
on a 2-D 5 nmx 5 nmx thickness grid using an exchange 
stiffness constant of 1.05 x 10 _11 J/m, with Ms values be- 
tween 700kA/m and 800k A/m and either 20 nm or 40 nm 
thickness. All calculations with the model used an ex- 
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change length of 5.85 nm and Ms values matching the 
simulations. 

Comparison of the m(h) = M(H/M S ) /(n M s V) and 
r(h) = Ar(H/Ms) /R curves are shown in Figures[2]a and 
b. Clearly the susceptibility-corrected 3rd order RVM 
provides the best estimate of both magnetization and 
vortex displacement as a function of field, while the un- 
corrected version exhibits the poorest performance. Both 
the DVPM and TVM provide good estimates of vortex 
position with field for displacements up to i?/2, but only 
the DVPM simultaneously gives a good description of the 
magnetization. 

Two other metrics have been applied to evaluate the 
performance of the analytical models near zero field in 
past work: initial susceptibility, and the frequency of the 
lowest order excitation mode of gyrotropic vortex motion. 
Both of these parameters primarily depend on the aspect 
ratio of the disks. Initial susceptibility is easily calculated 
from m (h) for each model. Using the collective coordi- 
nate approach^, it may be shown that the gyrotropic 
mode frequency is f Q = k/2i:G where G = 2ttLMs /jo 
with 7 = 1.76 x lO 11 ^ 1 !^ 1 , and dE tot /dr = kJ^. For 
the RVM and TVM, k is f3 and 4a respectively. For the 
piecewise combined model, k may be computed in the 
unpinned and zero field case 

_ 4(q + 7C 2 /?)/?a 

Dynamic simulations were performed using a realistic 
damping factor (0.02) but otherwise matched the param- 
eters used in the previous simulations. The poor perfor- 
mances for magnetization description of the TVM and 
uncorrected RVM manifest as incorrect estimations of 
the initial susceptibility. However, both approach the 
simulation results for squat disks, corroborating previous 
results^^ and demonstrating the general utility of these 
models. By comparison, the DVPM and corrected RVM 
provide excellent estimates for all aspect ratios investi- 
gated. Previously, only the TVM has provided reason- 
able estimates for the gyrotropic frequency of the vortex 
state while the RVM has provided poor estimates. The 
success of the TVM is reproduced here, as is the failure of 
the uncorrected RVM. The susceptibility correction im- 
proves the RVM prediction, however it fails to match the 
performance of the TVM in the prediction of f Q . How- 
ever, the DVPM provides comparable performance to the 
TVM for low aspect ratios and improved performance 
with more squat ideal disks. 

IV. APPLICATION TO PINNING 

Having demonstrated the performance of the piecewise 
model in a perfect disk, pinning may now be considered. 
Adding pinning to the models is accomplished by adding 
functions of the form E p (b + R\a/2R — X p ) for a pinning 
site located at X p to equation^} or of the form E p (b—Xp) 
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FIG. 3. (Color available in online version.) (a) The 2-D 
energy of the DVPM (color scale) is plotted as a function pa- 
rameters b' (rigid) and a (flexible) for a small applied field. A 
single pinning site is included, and the pinning site coordinate 
axes i and j are shown. At the applied field, local minima ex- 
ist outside the pinning site (green disk at bottom) and inside 
the pinning site (blue disk at top) separated by an energy bar- 
rier lowest at a saddle point between the two minima (red disk 
in middle). The solid lines indicate the field evolution of the 
pinned and unpinned minima. Note that flexing evolves non- 
monotonically. (b) Solving for the minima as a function of 
field, and computing the corresponding magnetization yields 
an m-h curve that includes pinning effects (green and blue 
solid lines). The simulated results (black line with dots) for 
a 1 /im diameter, 40 nm thick disk with a single large pinning 
site at 105 nm from the center. The pinning site is included 
in the model calculation at 107.5 nm with an energy profile 
estimated from the simulation (inset). The model m-h curves 
extend over the field range where a local minimum exists, (c) 
Incorporating pinning into the RVM (dotted blue and green 
lines) results in much larger positional error (> 10 nm, shown 
uncorrected here) and a large underestimate of the hysteresis 
loop width. Inset is an image of the Ms suppression included 
in the simulation. The width of the pinning site is 40 nm, and 
at the center Ms =550kA/m. (d) Plotting the simulated, 
and computed DVPM and RVM magnetization as a function 
of displacement shows that the DVPM accurately captures 
the change in the character of the displacement of the vortex 
during pinning as well as the contribution of flexing to the 
magnetization. 



for the RVM in equation [T] For the RVM case, simply 
solving for the minima in energy permits a full solution 
of the problem. For the DVPM, the 2-D optimization 
required makes the problem more complicated, however 
this is critical to the success of the model. Plotting the 
pinning energy, pinning sites appear as linear troughs in 
a-b space (Figure [3] a). This permits a simplification of 
the optimization process by the consideration of pinning 
site coordinates defined by b — isin{6) + j cos(0) and 
a = icos(0) — j sin(0) for = tan~ 1 (2R/ R±). Switching 
to i and j coordinates allows independent minimization 
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and simplifies the problem. Computation of the values 
of b and a to minimize the function permits description 
of the quasistatic pinned behaviour, while locating the 
saddle points separating minima allows computation of 
the energy barriers separating bistable states. 

The critical feature of this minimization process is 
that the coordinate a may evolve non-monotonically 
with increasing field (Supplementary Material Movie M3) 
matching the qualitative non-monotonic evolution of the 
flexing of the magnetization distribution visible in simu- 
lation (Supplementary Movies Ml and M2). 
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FIG. 4. (Color available in online version.) (a) The DVPM 
m-h result is compared to a simulation of a 500 nm diameter, 
40 nm thick disk with Ms =700kA/m and a single pinning 
site located 80 nm away from the center, (b) The computed 
displacements are compared for the same simulation. For this 
comparison Ri =80 nm provides the best estimate of hystere- 
sis loop width. The inset is a plot of the optimal Ri value 
found by comparison to simulation as a function of aspect ra- 
tio. The blue fit line is included as a guide to the eye. (c) 
Comparison to a simulated 1500 nm diameter, 40 nm thick 
disk with a single pinning site at 200 nm from center shows 
that the DVPM begins to underestimate the magnetization 
as disk size increases, (d) The computation of the vortex po- 
sition and hysteresis loop width remains accurate. A Ri value 
of 375 nm was used. 

Micromagnetic simulations were used evaluate the pin- 
ning performance of the DVPM and, for comparison, the 
corrected RVM. The same simulation parameters were 
used as in the previous simulations. Pinning sites are 
mimicked using approximately circular regions of de- 
pressed saturation magnetization to modify the energy 
landscape of the disk (Figure[3]b inset). This leads to two 
contributions to pinning energy, the reduced exchange 
energy of the core in the low Ms region, as well as re- 
duced demagnetization energy when the core is centered 
on the site. The energetic profile of the pinning site can 
be approximated by considering the convolution of a 2- 
D Gaussian at various offsets with the profile of the Ms 
variation (Fig pp, inset). The Gaussian effectively ap- 



proximates the exchange energy density of the core, as 
well as the M z profile, providing an estimate of how the 
two energy contributions change as the core shifts rela- 
tive to the pinning site. Here a full width half max of 
17.2 nm is used for the exchange energy density. 

The performance can be evaluated by three metrics, 
the pinning site position error, the width of the minor 
hysteresis loops associated with pinning and depinning, 
and the combined computed pinned differential magnetic 
and positional susceptibilities. Figure 3 b-d show results 
for a 1 micron diameter 40 nm thick disk compared to 
the DVPM and the 3rd order RVM. The DVPM accu- 
rately captures both differential susceptibilities while the 
RVM fails to capture the positional slope. Both models 
feature effective position shifts of the pinning site. The 
DVPM agrees best with the simulation for a pinning site 
shifted 2.5 nm further from center than the actual sim- 
ulation, while a shift greater than 10 nm is best for the 
RVM. Most importantly, the computed entrance and exit 
hysteresis loops agree closely for the DVPM, but are al- 
most non-existent for the RVM. The deformability of the 
DVPM permits the vortex to flex ahead into the site, 
and linger in the site at a lower energy cost than the 
rigid model. 

Disk sizes between 500 nm diameter/40 nm thick and 
2000 nm diameter/20 nm thick were simulated with iden- 
tical Ms variation pinning sites. Unlike the perfect disk 
case, model performance for pinning is significantly af- 
fected by the choice of R\. The amount of flexing per- 
mitted affects the vortex position and magnetization in 
the pinning site, and consequently, the energy barriers 
for entrance and exit from the pinning sites. The op- 
timal R\ reflects the overall rigidity of the disk behav- 
ior. In general the value R\ = R/2 provides reasonable 
results. However for disks below 1/im in diameter (for 
40 nm thickness) reduced R\ values provide better pin- 
ning performance (Figure[3]i, inset) reflecting the increas- 
ing rigidity of smaller disks. The DVPM was found to 
give good estimates of hysteresis width and vortex posi- 
tion for all disk sizes when an optimized value of R\ was 
used. For disks significantly larger than 1 /am in diame- 
ter, at 40 nm thickness, the pinned magnetic differential 
susceptibility was found to be underestimated, reflecting 
a change in the character of flexing in larger disk sizes. 
Figure 4 shows a 500 nm and 1500 nm diameter result for 
comparison. 



A. Two Dimensional Pinning Potentials 

The DVPM provides excellent performance in the de- 
scription of ideal disk behavior and pinning for idealized 
simulations. However, in application to real samples, the 
treatment of pinning sites located directly along the path- 
way followed by the vortex in the absence of pinning as it 
is deflected by field is limiting. As noted in recent numer- 
ical simulation work on pinning 39 , a more realistic case 
is to consider pinning sites near, but not centered on, the 
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Reduced Field (H/M s ) Ax Displacement (nm) 

FIG. 5. (Color available in online version.) (a) A 3-D plot 
shows the combined pinning site potential and harmonic po- 
tential for deviations orthogonal to the path defined by the ap- 
plied field. The computed path of minimum energy, Ax(Ar) 
is plotted as a red thread, (b) At top the 2-D pinning poten- 
tial is plotted with the computed minimum energy path (red 
line). Below, the equivalent 1-D potential is presented. The 
equivalent potential incorporates contributions from both the 
pinning potential and harmonic trough. The color gradient 
on the line matches the color scale in panel (a), (c) The 
magnetization curve computed from the potential in (b) is 
compared against a simulated curve incorporating Ms sup- 
pressed regions with the same 2-D distribution. The depth of 
the 2-D potential used in the model calculations is estimated 
from the simulation. Agreement is close, though some devia- 
tions show up as the vortex displacement increases and a large 
energy change is encountered, (d) The simulated Ax(Ar) is 
compared to the computed minimum path showing excellent 
agreement. 

field defined path. This can be incorporated into the 1-D 
model presented here by computing the 1-D equivalent 
potential of the actual 2-D path followed by the vortex. 
Deviations orthogonal to the path defined by the applied 
field have an energy cost approximated by nAx 2 where 
the value of k is given by equation [4j Since the mag- 
netization induced by these deviations is orthogonal to 
the applied field, the energy is effectively static and can 
be summed with a 2-D distribution of pinning sites to 
form a trough guiding the vortex through the 2-D energy 
landscape. It is then possible to compute the minimum 
energy pathway Ax Q (Ar) that the vortex will follow as it 
is deflected (Figure pk). Computing the total static en- 
ergy, pinning plus the trough energy, E(Ax Q ) yields an 



equivalent 1-D potential as a function of Ar (Figure^). 
This potential can then be summed, as the Gaussian pin- 
ning sites were previously, with the potential for a perfect 
disk including field. Solving for minima as before allows 
computation of the evolution of the magnetization and 
vortex position in the 2-D potential. 

This approach is applied to a simulation that incorpo- 
rates a 2-D distribution of 10 nm diameter pinning sites 
with various values of suppressed Ms near the field de- 
fined path. As before, the pinning sites are incorporated 
into the model as Gaussian wells with depths estimated 
from the simulation and profiles computed by convolving 
a Gaussian with the profile of the Ms variation. The 
computed 2-D path agrees well with the vortex position 
extracted from simulation, as does the computed magne- 
tization (Figure [2] c and d). Some disagreement is noted 
as the deflection increases close to the effective R/2 limit 
of the model, and the vortex passes over a large bar- 
rier. In this computation, a sparse distribution of 2-D 
sites ensures a unique Ax G (Ar), however, in principle, 
this approach can be extended to bistable states in Ax 
by consideration of multiple vortex tracks. Computation 
of the energy barriers separating tracks, however, would 
require a more complete minimization. 



V. CONCLUSION 

The piecewise approach applied to develop the DVPM 
yields a highly functional analytic model that makes 
quantitatively accurate predictions of a wide variety of 
properties of a vortex in a disk. Most notably it pro- 
vides a powerful description of vortex pinning and pro- 
vides greater physical insight into the behavior of the 
vortex during pinning. The model holds promise as a 
tool in probing the modification of pinning in techno- 
logically pertinent thin films to better understand effects 
such as ion damage, while the piecewise approach demon- 
strated may, in future, be generalized to other geometries, 
permitting quantitative computation of device behavior 
without cumbersome simulation. 
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VI. SUPPLEMENTARY MATERIAL: A 
DEFORMABLE MODEL FOR MAGNETIC 
VORTEX PINNING 

The first movie (Ml) shows the magnetization distri- 
bution of the disk from the simulation that generated the 
data used in Figure 3 in the main text . The magneti- 
zation direction is indicated by the color. Red indicates 
magnetization in the positive x direction, green is nega- 
tive blue is positive y and yellow is negative y. 

As the vortex is displaced by the field, the magneti- 
zation distribution warps away from the circularly sym- 
metric initial vortex state. As the vortex interacts with 
the pinning site, the warping changes, and the flexing 
of the distribution decreases. This results from the con- 
tinued biasing of the magnetization distribution outside 
of the core effectively mimicking continued displacement 
of the outer region while the core has reduced mobility. 
This can be challenging to see on the color scale. A sec- 
ond movie (M2) of the same simulation shows only the x 
component of the magnetization in a contour plot. The 
non-monotonic evolution of the contours stands out more 
clearly. 



The third movie (M3) shows the evolution of the en- 
ergy (color scale) of the vortex for a given location depen- 
dent on the parameters a' and b' computed in the model 
as magnetic field increases. The vertical axis is a! and b' is 
the horizontal coordinate. The single pinning site shows 
up as a trough across a' \b' space. As the field increases, 
the unpinned vortex position (blue dot) increases and 
a second minimum appears in the pinning trough (red 
dot). As the field increases further, the vortex position 
outside the pinning site continues to increase in both ol 
and b' but the pinned minimum decreases in a'. This 
effectively captures the behavior in the simulated movie. 
The flexing decreases, while the outer region continues to 
displace. 



In addition to providing the opportunity to determine 
the positions of energetic minima to compute the un- 
pinned and pinned magnetization and vortex displace- 
ment curve as a function of field, the saddle point sepa- 
rating the two minima may also be computed, allowing 
calculation of the energetic barrier separating the two. 



